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Abstract 

We make an attempt to discuss in detail the effects originating from the 
final state interaction in the processes involving production of unstable ele- 
mentary particles and their subsequent decay. Two complementary scenarios 
are considered: the single resonance production and the production of two 
resonances. We argue that part of the corrections due to the final state inter- 
action can be connected with the Coulomb phases of the involved charge par- 
ticles; the presence of the unstable particle in the problem makes the Coulomb 
phase "visible" . It is shown how corrections due to the final state interaction 
disappear when one proceeds to the total cross-sections. We derive one-loop 
non-factorizable radiative corrections to the lowest order matrix element of 
both single and double resonance production. We discuss how the infrared 
limit of the theories with the unstable particles is modified. In conclusion we 
briefly discuss our results in the context of the forthcoming experiments on 
the W+W- and the ti production at LEP 2 and NLC. 
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1 Introduction 



Processes involving production of unstable fundamental particles such as top quark, 
electroweak massive gauge bosons, Higgs boson are now at the frontier of both the- 
oretical and experimental high-energy physics. Investigation of such processes can 
provide a valuable information about fundamental parameters (mass, width) of the 
heavy unstable particles. Therefore the theory has to give reliable predictions to 
meet expected experimental precision. However it must be recognized that the ac- 
curacy of the theoretical description of the processes involving unstable particles and 
the planning accuracy of the measurements have been never previously combined. 
There exists a number of recent theoretical inventions which all are connected with 
the accurate description of the unstable particle in the vicinity of its pole: the S- 
matrix approach to the Z-pole [|l| and the extension of this scheme to a more general 
cases [0, colour rearrangement phenomena [Q, etc. 

Currently, the main source of our experience in the field is the Z-boson physics. 
However, as it has been mentioned in and will be quite clear from the forthcoming 
discussion, Z-pole description is distinguished due to the following points: 

1. Z-boson is neutral; 

2. Main results have been obtained for the process e'^e'^ ^ Z ff, i.e. for the 
production of the Z-boson in the s-channel. 

These features greatly simplify precise description of the Z-boson production cross- 
section. 

The basics of the theoretical approach to the processes with the unstable particles 
can be described as follows: when two unstable particles are produced all Feynman 
graphs can be divided in two classes: the first one includes the graphs without 
interactions between decay products of different unstable particles, while the second 
includes such graphs in which two decay processes are not independent. Generally 
speaking the graphs of the second type give a correction of the relative order aT/Mi, 
where a is an appropriate (depending on the process) coupling constant and Mi is 
the characteristic scale for the momentum flow inside the loop. For example the 
graphs with the Z-boson exchange between decay products of unstable particles 
provide Mi ^ Mz- Hence their contribution is negligible. On the other hand for the 
graphs with photon or gluon exchange the scale Mi is of the order of the width of 
the unstable particle F. Consequently such graphs have no additional suppression 
in comparison with the factorizable ones []. It is not difficult to convince oneself 
that the dominant contribution comes from the soft photon or gluon region [4-6]. 
This contribution is a non-trivial function of the width of the unstable particle. The 
reason is that the soft massless particles probe the mass-shell limit of the theory and 
it is the width of the particle that changes the mass-shell behaviour of the resonance 
in comparison with the stable particle. Hence in this problem (similar to the case 
of the threshold production of the unstable particles 0]) one have to use resummed 
propagator of the unstable particle. 

^As for gluon exchange, this also verifies the use of the perturbative QCD for the calculation 
of these contributions in the reactions with the top quark(s), since Mioop ~ sets the scale on 
which as must be evaluated for these non-factorizable contributions. 
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The other connected physical problem which has attracted serious attention in 
the past years is the photon or gluon radiation of the unstable particles 0. It 
turns out that in the specific kinematical configuration the radiation is completely 
insensitive to the decay. The reason for this can be found in the conservation of the 
charged or colored currents. 

As is well known the soft physics generally "suffers" from the cancellation be- 
tween real and virtual corrections. It is therefore quite desirable to clarify how 
and when this cancellation occurs when unstable particles are involved and to what 
extent one can study real and virtual corrections independently. 

Recently some progress has been achieved in the understanding of these prob- 
lems. It has been argued [4-6] that non-factorizable corrections do not contribute to 
inclusive quantities, e.g. to the total cross-sections, when both virtual and real ones 
are taken into account. This can be viewed as the extension of the Bloch-Nordsieck- 
Lee- Nauenberg-Kinoshita cancellation for the processes involving unstable particles 
0.0 However our analyses shows that it is not only usual real-virtual cancellation 
but something more involved. This fact becomes clear if one studies the influence of 
the non-factorizable radiative corrections on differential distributions. These distri- 
butions (say in the invariant mass of the decay products of a resonance) are a possible 
tool to investigate fundamental parameters of the unstable particle (for recent dis- 
cussion of the ly-boson and the top physics see refs.0, [T^ and references there 
in). It turns out that differential distributions are affected by this non-factorizable 
interactions (to the best of our knowledge this fact has been flrst noted in the ref. 



11 1 for the speciflc case of the top threshold production). 

In what follows we evaluate 0(a, a.^) non-factorizable resonance radiative correc- 
tions to the differential distributions in the invariant mass of the unstable particle(s). 
If the integration over invariant masses is performed, these radiative corrections dis- 
appear 0,0. Calculating radiative corrections to the differential cross-section we 
are able to clarify the physical origin of this "inclusive zero" . 

However, it turns out that the shape of the differential distributions and the 
position of the maximum of the differential distribution in the invariant mass of 
the resonance which could be naively identifled with the pole mass of the unstable 
particle are affected in the energy region slightly above threshold of two resonances. 
For the particular case of two W^-bosons this region is approximately 170 — 190 GeV 
being almost the same as the LEP2 energy region. 

Absolutely deliberately we do not consider the real threshold region (i.e ^/s — 
2M ~ 0(r)). In this region new physical phenomena appear (bound state forma- 
tion, etc.) and our analyses would be more complicated there. Our idea is to get 
the most clean laboratory for the effects which are completely connected with the 
unstable nature of an appropriate particle. Threshold region represents a special 
case and has to be discussed separately (for the top threshold production see |]rT|, 

0)- 

Let us note that through out the paper we use the Breit-Wigner with a constant 



^ Note that in general this cancellation differs from the cancellation known from the Z-pole 
physics. This fact clearly follows from the ref. [Q. Indeed, in describing Z-pole we deal only with 
the initial-final non-factorizable interaction which is much more simple. In general case there is 
also final-final state interaction which brings some new features to the problem. 

■^We remind that the lowest order differential distribution is the common Breit-Wigner. 
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width as the propagator for the unstable particle. As we are concerned with the 
corrections of the order 0{a) to the lowest order result we can safely use the lowest 
order propagators for the unstable particles since all modifications show up only in 
higher orders. 

Subsequent part of the paper is organised as follows: next section is devoted to 
the investigation of the single resonance production; more involved scenario with the 
production of two resonances is discussed in the section 3 where all basic formulae 
are presented. In the section 4 we analyse our results in a more informal way. 
Conclusion of the whole work is given in the section 5. A number of helpful formulae 
are presented in the Appendix. 



2 Simple model 

For simplicity we start with the model describing scalar particles which interact 
with the "photon" field. Suppose one of this particles (we call it W) can decay to 
two other (electron and neutrino for simplicity). Our W particle is produced by 
some neutral current (virtual photon) together with the other stable particle [B- 
particle). In such a model W and B have opposite electric charges. 

In what follows we discuss reaction 7* —>■ W^B~ e^vB~ taking into account 
0{a) non-factorizable radiative correction^. 

We consider this process in the center of mass frame of the virtual photon. Then 
it carries the total energy ^/s and the zero three momentum. We are interested in 
the distribution over the invariant mass of the W particle. In order to describe this 
distribution we introduce a parameter: 

where My/ is the pole mass of the W and is the invariant mass of the final eu 
system. The Born graph is shown in the fig.l. Above the production threshold of 
the W particle the Born graph has the resonant propagator forcing produced W to 
be almost on shell. Non-factorizable virtual corrections are also shown in the fig.l. 
Let us first discuss the graph with the B e interaction. Since we are interested in 
the corrections of the order 0(a) we have to get the resonance denominator from 
this graph. Consequently the loop momentum must be small in order not to shift 
the VT-particle propagator far from the pole. From this it is clear that the only 
loop momentum region which can provide such "resonance" correction is the soft 
region, where one can use soft-photon approximation (cf. ref. 0]). In the soft 
photon approximation the amplitude of this process reads (we use the Feynman 
gauge through out the paper): 

ATxiaMQ f d^k Ape ■ Pe _ , ^ 

D{py,) J {2n)Hk' + ^e){2psk + ^e){2p,k -^"^^^^ 



It 



Here Mq is the Born amplitude and 



^'^Pw) = — ^ I (3) 



"^In some sense this case corresponds to the process t — > W~^b. 
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is the propagator of the W particle with the finite width inchidcd exphcitly. 

To perform the integration over k we first integrate over zero component of this 
four vector. There are four poles in the complex plane of the ko variable. Let us 
integrate over the lower half plane of the A;o~complex plane. Then two poles have 
to be taken into account: one from the B particle propagator (" the particle pole") 
and the other one from the virtual photon propagator (" the photon pole"). It can 
be seen that in contrast to the soft photon approximation in the QED with stable 
particles, in our case the contribution due to the virtual photon pole does not cancel 
corresponding soft photon emission immediately. However, their difference appears 
to be pure imaginary and hence does not infiuence differential distributions (for 
more detailed discussion see section 4.2). Keeping this in mind we take into account 
the -B-particle pole only. 

The contribution due to the B particle pole reads: 

M„ = -Ana ^iRlllJp- /, (4) 
AE.Eb D{pw) ^ ' 



~ I T^T( — — T^wl \ ^ ^^^"^ ~ ^)Uo=kvB- 

J (27r) ((Vfik) -\s?) (kvB-kVe-^e) 



Here the quantity Vj is the three velocity of the i-th particle. Using momentum 
conservation we get for the W propagator: 

D{pw - fc)|fco=kvB = ^- 7=^= -• (5) 

Mw{5w + i^w - 2^s/M^kvB) 

In order to compute residual integral over k it is useful to exponentiate the propa- 
gators, introducing two different "times" . The amphtude reads then: 

-A-naApB ■ PeMo r _^ -tw f exp{ikr(t, r)} 

AE,EbMwD{pw) J ^ (27r)' (vBk)'-k2 

where r{t, r) stands for: 



r(t,r) = (ve-v^)t-2^^v^r. 

The integral over k is recognized to be retarded Coulomb potential of the particle 
moving with the velocity and hence the result of k integration can be found in 
the text-books on classical electrodynamics: 

, , , , /■ d^k expjikrj 1 

Hr) = -47r / { = 6 

Here ns is the unit vector parallel to the B-particle velocity and rj_ is the component 
of the vector r transverse to the vector n^. The final expression which can be 
obtained in this way is: 
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Let us note that the way we proceed is quite similar to the eikonal approximation 
for the high-energy scattering. It is well known in that case and can be proved in 
ours that the leading contributions from the eikonal graphs to the amplitude can 



be summed up . The result is the Coulomb phase [|T4| of the wave function of 
the charged particle. Usually the Coulomb phase is not important due to its pure 
imaginary nature. We shall see that in our example this is not the case and that the 
residual contribution from the Coulomb phase survives in the final result. 

Integrating the last equation over t and r we neglect the terms which are pure 
imaginary and hence do not contribute to the differential cross-section at the 0{a) 
order. We get then: 

Ma = -a I log(- — 1. (8) 

In this equation Ve_L is the component of the vector Ve transverse to the vector n^. 
Note that the factor in front of the logarithm is nothing but the Lorentz boosted 
Coulomb factor a/|vi — V2|. This factor has the following limits: when velocities 
are small it turns to a/jv^ — ^e\ hence reproducing usual expansion parameter for 
the Coulomb problem |T3 while in the limit |ve| — >■ 1 or |vb| — 1 this factor equals 
to a and hence appears to be independent from the kinematic of the process. 

Let us now discuss the photon exchange between W and B (see fig.l). On the first 
glance this graph does not seem to be non-factorizable correction we are interested 
in. However gauge invariance arguments do not allow us to exclude this graph from 
the consideration. We study this graph in the soft-photon approximation neglecting 
the contribution of the photon pole (see sect. 3. 2). Calculation is quite similar to the 
previous one and results in the following contribution to the amplitude: 

Mb = a Mo t log I — . (9) 

There are no other corrections of the non-factorizable origin which influence differ- 
ential distributions. For example the interaction of the W with the electron is of 
the initial-final state interaction type and hence has rather simple pole structure. 
The infrared contribution from this graph is completely cancelled by the corre- 
sponding real emission ( see also discussion in the ref. []T3[ for the stable particle 



case). Hence all radiative corrections which are of the non-factorizable nature and 
are not cancelled by the emission of the soft photons are given by the sum of two 
amplitudes presented above. The sum of this amplitudes gives us the result for the 
non-factorizable corrections: 

Mn/fact = aMoi log ( J ( I 1 - , ^ )■ 10 



Let us write corresponding contribution to the cross-section in the following form: 

da = dao K = —2 rj arctg(:i; — ), (11) 

/ 1 - VwVb 1 - VeVs N 

f] = a { r . = . 



e± 
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Here dao is the lowest order cross section. The important point to be noted here is 
that in the relativistic hmit for this equation the cancellation between contributions 
from WB and Be interaction occurs. This resuh recovers the "non-observabihty" 
of the Coulomb phase. We could expect this compensation because in the ultra- 
relativistic limit the spectator ("B" -particle) does not distinguish transverse move- 
ment of the electron and hence does not notice that the charge movement has been 
changed. As the result, the Coulomb phases of the resonance and its decay products 
add coherently to a pure imaginary quantity reproducing the non-observability of 
the Coulomb phase. 

Following above discussion we can exponentiate our results given in the eqs. 
(12)-(13) (cf. ref. 0) to get: 

M^/f,,, = MoTil + tr])e''^o( \' ^^2) 

where is an imaginary phase which is not relevant. It is straightforward to find 
the contribution of this amplitude to the differential cross-section: 

da = dao —n — v exp{-2 77 arctg(-^)|. (13) 
sh['7ir]) Tw 

To see how this correction disappears when we proceed to the total cross section 
[4-6] we integrate previous equation over the range of the W masses (of course with 
the usual approximations in mind) and get 

On the one-loop level this cancellation is clearly visible since the a.Tctg{Sw/^w) is 
the odd function of the resonance off-shellness, while the usual Breit-Wigner is the 
even one. Therefore their convolution is zero. 

It is clear however that this correction influences differential distributions in the 
invariant mass of the produced resonance. To get an idea of what one gets in the 
realistic situation let us imagine that we deal with the production of two equal mass 
resonances and the integration over invariant mass of one of them has been already 
performed. We treat the resonance which is "integrated out" as a stable particle. 
The velocity of the "electron" is taken equal to unity. Then as it has been noted 
above the factor 77 turns out to be independent from the scattering angles of the 
final particles. In this case this factor reads: 

a (15) 



2/3 
where 



is the on-shell velocity of the produced resonance (as far as we are not too close to 
the threshold we can use this on-shell value for the velocity). We note here that 
the T^-factor in the eq.(|T3p very quickly goes to zero if the total energy increases. 
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For example, taking the mass of the resonance equal to 80 GeV one sees that for 
^/s = 170 GeV we have r] = 0.65 a and for y/s = 200 GeV this factor decreases up 
to 0.13 a suppressing this correction roughly to one order in magnitude. 

However, this correction influences differential distributions in the invariant mass 
of the ly-particle decay products moving the peak to the lower values of the reso- 
nance masses. The result for the corrected distribution is shown in the fig. 2 for the 
usual values of the width and the mass of the ly-boson and for different energies 
of the process. The position of the maximum of the distribution differs from the 
same quantity defined by the Breit-Wigner propagator. The position of the new 
maximum is: 



Though the pole position is not affected too much for the realistic values of the 
particle width and the coupling constant it is still comparable with the planned 
accuracy of the W mass determination in the intermediate energy region -y/i ~ 
170 — 190 GeV. For higher energies this corrections are strongly suppressed hence 
having no importance from the experimental point of view. 
Let us make some comments now. 

There exists the S-matrix approach for the description of the gauge boson pole 
which was originally proposed for the description of the s-channel production of the 
Z-boson ||l[. The basic idea of this method is to start from the analytical properties 
of a given amplitude. Our analysis shows that in order to apply this method to 
the charge boson production one must claim that there is a branching point in the 
complex plane of the invariant mass of the resonance but not a pole. Corresponding 
intercept of the branching point appears to be non-trivial function of the kinematic 
of the process. From this we conclude that a theoretical analysis of this situation 
will be more complicated and there is no straightforward extension of the S-matrix 
pole scheme to the processes involving charge unstable gauge boson(s) production^. 

Our next remark concerns the top decay width. As it is clear from the exact 
expression for the rj factor the non-factorizable corrections to the differential distri- 
butions over invariant mass of the ly-boson decay products are negligible due to 
the small mass of the 6-quark. 

Let us also outline how the calculation of the radiative corrections to the single 
resonance production must be performed. We stress once more that the radiative 
correction which is presented in the eq.(13) is the only one which is usually referred 
as non-factorizable. As we have also traced the cancellation of the real emission 
against virtual photon poles we can formulate the practical recipe for the calculation 
of radiative corrections to the differential cross-section of the reactions with one 
unstable particle: 

1) The leading order differential distribution is given by the known formula: 



where cxo is the on-shell cross section for the production of the particles W and B, 
and dr{W —>■ ev) is the on-shell differential partial width of the H^-particle. 
2) In order to compute 0{a) corrections to this formula one needs: 

^This fact and a possible modification of the pole scheme are discussed in the ref. 




(16) 



da{pl,) ^ ao(M^) Mw dVjW ^ ez/) 



(17) 
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• to substitute one-loop results for all quantities in the previous formula; 

• to add our result for the "Coulomb phase" contribution (eq.l3) to the above 
formula, since this is the only contribution of the order 0{a) which comes 
from the non-factorizable interaction. This prescription already takes into 
account partial cancellation of the non-factorizable real corrections against 
corresponding real ones. 

There is one subtle point in the preceeding discussion. One can get an idea 
that we make a double counting, i.e. we include the soft region of the triangle 
{bW interaction in the terms of the model) to the "narrow width approximation", 
while this region is also accounted in the calculation of the non-factorizable radiative 
corrections, which according to the recipe are added later by hands. We note in this 
respect that the soft region for the "on-shell" triangle is completely cancelled by the 
soft emission. In this case we do not get any contribution from this region because 
the Coulomb phase for the stable particles is pure imaginary and hence disappears 
from the observables. Consequently there is no double counting and our recipe is 
simple and reasonable. 



3 Production of two resonances 

Now we are in position to discuss similar problem for the case when two unstable 
particles are produced in an appropriate reaction. We note that non-factorizable 
radiative corrections to the processes involving production of two resonances have 
been discussed in the literature. Namely total cross- sect ions [4-6] and various distri- 



butions in the non-relativistic ( threshold) limit [|TT| have been analysed. As we have 
mentioned in the introduction, we do not discuss threshold region in what follows. 

Below we calculate double resonance 0{a) non-factorizable radiative corrections 
to the lowest order matrix element. Through out the paper we consistently neglect 
single-resonance and background contributions. 

Next important remark concerns logarithmic and polylogarithmic functions which 
appear in the result of this calculation. Generally we need to evaluate this functions 
in the complex plane. Hence it is important to fix conventions for the cuts of this 
functions. All logarithmic and polylogarithmic functions which appear in our final 
formulae have the usual cuts i.e for the logarithms it goes from to — oo along the 
real axis and for the polylogarithms from 1 to +cxo along the real axis. 

For concreteness we consider the process 7* — *■ — *■ bW^bW~ as the basis for 
further discussion. However, for the energy region sufficiently far from the threshold 
the results of our calculation appear to be general and are not restricted to a concrete 
process. 

It is worth to note from the very beginning that the results which one obtains 
for two resonances appear to be more complicated and are not so transparent from 
the physical point of view as compared to the case of a single resonance production. 
If it is possible we try to appeal to the physical picture rather than huge formulae. 
The generic graphs for the non-factorizable radiative corrections are presented in 
the fig.3. 
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3.1 Three-point function 

We start our consideration with the usual three-point function. As it is clear from the 
previous section the contribution of this graph is unavoidable due to the gauge invari- 
ance arguments. Corresponding amplitude diverges logarithmically in the soft-gluon 
approximation. To avoid this divergence we introduce the cut-off A and restrict the 
integration region to the values of the loop momenta < A^. This regularization 
is not Lorentz- invariant, so we use the center of mass frame everywhere. The result 
of the calculation is Lorentz invariant anyhow. 

As we work in the soft gluon approximation, we are interested in the contribution 
from the region /c ~ F. The natural requirement for the cut-off A is then 

r << A « y/s,Mw. 

When the momentum of the virtual or real gluon is much large than F we can 
neglect the width and the off-shellness in our formulae and work with the usual 
expression for the radiative corrections. The important point is that the usual 
radiative corrections (both real and virtual) drop out from the observables in the 
soft gluon approximation. Hence we expect that the cut-off A will not enter our 
final formulae. 

We mention here that the complete result for the three point function with 



two unstable particles is known in the closed form ||15[. One can use this result 
without any approximations avoiding the questions associated with the ultraviolet 
divergence of the soft gluon approximation. However, the soft part of this triangle 
can be correctly reproduced by the soft gluon approximation. This part is necessary 
for our analyses providing the complementary gauge invariant part for the four and 
five point functions. 

3.1.1 Particle poles — The amplitude for the process under discussion corrected 
due to the three-point function in the soft-gluon approximation is: 



i^k 1 



(27r)' {Di - 2pik) {D2 + 2p2k) 
Here 

Di=pi^-mt^ + imtTt, 2 = 1,2 (19) 

is the off-shellness parameter which we use further. 
Cp stands for the usual colour factor: 



A^c'-l 

2Nr ' 



Again, there are two "particle" and two "gluon" poles in the eg . (P^Sf ) . We perform 
the integration over the lower half of the complex plane. Let us discuss the particle 
pole first. 

Taking the particle residue in the ko complex plane, introducing the cut-off 
and splitting momentum integration into parallel and perpendicular components 
with respect to the resonance velocity ( in the center of mass frame two produced 
resonances move in the opposite directions) we get 

M,,rt = Mo^{l + P') I, 

TV 
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(271) J (2'Kf {D2/ E + P ~ - kl iDi + D2)/E + P k. 

Here E = ^/s is the total energy of the process and [3 is the velocity of the particle 
with the energy E /2 and the mass mj. 

The integration over kz is restricted to the region —A < kz < A and the integra- 
tion over k± to the region k'^ < — fc^. 

It is seen from the eq.(22) that the integration over k± is logarithmic and hence 
can be performed immediately. Integrating by parts over kz we get: 



4/3 J ^ ^ Errit rrit 
2(3{D2/E + (3kz) 2(3{D2lE + (3 kz) ~ 2kz 



Mpart = -Mo — — dkzlogi—- + )P{,kz), (21) 

{D^/E + pkzf-k^ {D2/E + [3kzf-kl ■ 

The simplest way to evaluate this integral is to use analytical properties of the 
integrand in the complex plane of kz variable. The logarithmic function has branch- 
ing point below the integration path and the singularities of the function P{kz) are 
simple poles. Using Cauchy's theorem we rewrite eg. (pi)) as an integral over the 
half-circle of the radius A in the upper complex half-plane, taking into account the 
residues where necessary. The result of the integration is then|: 

M,art = -Mo^^-^^p-{2n I log(e(-l)) + vr^) (22) 

where the function S,{x) will be used further through out the paper. This function 
reads explicitly: 

i{x) = (1 + /3a;) ^ + (1 - [3x) % (23) 

This result shows a peculiar property - after the choice of the particular inte- 
gration contour we have lost the symmetry between two resonances, in spite of the 
fact that the original integral has such a symmetry. As we have learned before, the 
particle pole gives the Coulomb phase, hence the absence of the symmetry in the 
particle pole contribution means that a part of the Coulomb phase is hidden in the 
contribution due to gluon pole. 

Let us transform the integral to the coordinate space to study the space-time 
picture. Doing so we recognize that the eikonal approximation provides simple 
deterministic picture: a particle is moving with the constant velocity in the constant 
direction, the "density matrix" is the moving delta-function. In the case when the 
particle is unstable everything is the same, except the normalization of the density 
matrix - it is not a constant anymore. 

Starting from the expression for the amplitude presented in the eq. ([I8| ) we intro- 
duce Schwinger-Fock proper time for each of the resonances. The integration over 
loop momentum reduces to the evaluation of the Fourier transform of the gluon 
propagator to the coordinate space. The result is well-known: 



/ 



(27r)^P + 2e An'^x'^-ie' 



^We again neglect all pure imaginary quantities in this result. 
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In our case the four-coordinate of the gluon propagator is the difference in Lorentz 
coordinates of the resonances: 



In the center of mass frame we rewrite the result in the following way: 

M = -Mo— — (1 + /? ) / dTidr2- -2 — ■ 72 r- (25) 

^ J (T1-T2) -I3\ti+T2) -^e 

This expression exhibits poles on the integration path. The position of this poles 
corresponds to the movement of one particle in the field produced by the other when 
retardation effects are taken into account. The residues in these poles provide us 
with ( we again drop all pure imaginary quantities): 

M,„,e = Mo ^^i^(^7rlog(e(-l))+^7rlog(e(l)) + 7^2). (26) 

This expression has all desirable symmetry properties and corresponds to the Coulomb 
phases of two resonances which they acquire in the field of their partners. 

As the space-time picture shows that our understanding of the Coulomb effects 
is still valid we proceed further and extract the residual Coulomb phase contribution 
from the gluon pole. We will not use the proper time representation systematically 
and continue evaluation of the three point function in the momentum space. 

3.1.2 Gluon poles — So far wc have studied the "particle" pole contribution to 
the amplitude. Now we are in position to discuss the contribution of the gluon pole. 
The separation of particle and gluon poles in our calculation is useful due to the 
fact that the contribution from the gluon pole of the virtual graph is in very close 
analogy to the corresponding real emission. Hence if we get the contribution of the 
virtual gluon pole it is a matter of machinery substitutions to obtain the amplitude 
for the soft real emission. 

Taking the residue of the gluon propagator, which is located in the lower half- 
plane of the /cq variable, we find that the integration over transverse component of 
the loop momentum is again logarithmic and hence straightforward. We perform 
one integration by parts and arrive finally to the following representation for the 
gluon pole contribution: 

M, = Mo^^^^^iAr + A2 + AS + A,) (27) 

where 

/D1 + D2 , 2/5A;,w f3 (5 



J ^ hirrii rrii ■' ^ 



mt ^^Di/E - A + pk^ D2/E + A + l3k, 

f Di + D2 , 2pk, ^^ 1-/3 l±^__\ 

'-J/"' ''^^ Emt ^ mt )yD,/E-{l-P)k, D,/E + {1 + P)kJ' 



A 



A, = j dkAh) log + ^) U/E+(l + 3)k. + D2/E+(l + 3)kJ ' 



mt ^^Di/E+{l + (3)k, D2/E+{l + p)k, 
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A,= dkA-h) log [-——- + t;-) ( n _ n - + 



-A 



Emt mt ^^D^/E-{1-P)h D2/E-{l~p)h 



Let us discuss the advantages of this representation: evaluation of the integral Ai 
can be immediately reduced to the integration over the semi-circle of the radius 
A, as a consequence it will not depend on the off-shellness and the width of the 
resonances. Therefore it will be completely cancelled by the real emission. The A2 
term is the extracted contribution of the particle pole ( hidden Coulomb phase, as 
it has been called above), the last two terms are specific for the gluon pole. The 
calculation of this integrals is straightforward due to the fact that all of them are 
of a polylogarithmic type. It is clear that we need to evaluate polylogarithms and 
logarithms as the functions of the complex argument. We note in this respect that 
all important points for performing logarithmic and polylogarithmic integrals in the 



complex plane have been discussed long ago in the ref. |16 . 

As has already been mentioned, the contribution from the Ai term is completely 
canceled by the real emission hence we do not present it here. The result for the A2 
reads: 

A, = -n^-2m log(e(l)). (28) 

The most involved is the evaluation of both ^3 and A4 contributions. The result 
which one obtains after direct integration is: 

+ = l0g(C)l0g( ^,^ Tl^?? ) - ^log^(Ml2|) -l0g(Mi2|)l0g(^ 

\zi\\z2\mt ' 2 ^1^2! 



2Li2(C) - Li2(-Crfi2) - Li2( - + ^ 

^ "12 ^ 



+ ^(01 - 02)^ - (^^2 - '^l)(02 - 0l) - 7r(z/i + Z/2)- (29) 

Here the following notations are used: 

2;i=e(l), Z2=i{-\). 
Vi = arg(zi), (pi = arg(A), 

In the presentation of this result we split the answer into the modulus and the phase 
parts, and write each of them in a way which allows straightforward investigation 
of the /9 — > 1 limit. 

The result for the three-point function is then: 

Mu = Mpart + A2 + A^ + A4. (30) 

Here Mpart is defined in the eq. (|2^) . 

Now let us discuss how corresponding real emission can be obtained from these 
quantities. In particular we mean the interference of the gluons emitted by different 
resonances. It is straightforward to write the contribution of this interference term 
to the differential cross-section in the soft-gluon approximation. Direct examination 
of the momentum integral shows that it is sufficient to perform the following modi- 
fications in the result obtained for the virtual gluon pole to get a contribution due 
to the real emission: 
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• change the sign of the result. 

It is important to note here that this transformation does not influence analytical 
properties of the amplitude, hence we can perform it in the final result. It is seen 
from the eqs. (30-31) that the virtual pole contribution is not invariant under this 
transformation, hence the real emission will not cancel the contribution of the virtual 
gluon pole for the three-point function. 

To demonstrate this point we study the limit /? — > 1. It is straightforward to 
obtain the following from the eqs. (24), (30-32): 

• Particle pole: 27r02; 

• Gluon pole : 27r0i — 7r(0i + ^2) — \{<t>i ~ ^2)^ + const] 

• Virtual correction = Particle pole + Gluon pole. 

Transition to the real emission discussed above reduces to the transformation 0i — >■ 
TT — 01. As a result the sum of the real emission and the virtual correction in the 
limit /3 — > 1 equals to: 

7r(0i + ^2) + 20102 + const. (31) 

The constant term is independent from widths and off-shellnesses of the resonances 
and we do not present it here. The first term is the Coulomb phase of two resonances 
in the limit /3 — > 1 and the second is a correlation between the phases of two 
resonances. 

3.2 Four point function 

As a next step we consider the graphs with the gluon exchange between th or th. 
Evidently there is a symmetry between these two and having the result for one of 
them it is straightforward to reconstruct it for the other. For concreteness we study 
the interaction between t and h. 

3.2.1 Pcirticle poles — We start with the discussion of the particle poles. In 
this case it is not so easy to apply direct integration discussed in respect with the 
evaluation of the three point function and we use the following trick to reduce the 
necessary amount of work: in the soft approximation the product of two propagators 
of the unstable particles (cf. eq.(3)) can be decomposed as: 

Di^,-k)Dij., + k)^j^-^±^{^^^ + ^-l^). (32) 

Here p is the three momentum of the produced unstable particles ( p = Pi = — P2)- 
Examining the poles in the complex plane one finds that by appropriate choice of 
the integration contour the second term in this decomposition gives no particle pole 
contribution while for the first one it is sufficient to take the pole corresponding to 
the h propagator. In fact this decomposition leads to a mixture of the poles of the 
original expression. Hence strictly speaking the particle poles which are discussed 
below are some combinations of the original particle and gluon poles. However as we 
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have seen during discussing the triangle graph particle and gluon poles are hardly 
separated when we deal with the amplitude involving two unstable particles. So the 
"names" here are just a matter of taste. 

Taking the residue of the b propagator we get: 

ATrasCp 1 — V1V4 ^ 

Mpart = ^0— ^ ^ /, 

£k 1 

(27r)3 ((V4k)2 - k2)((Di + D2)/E + 2vik)(Di/E - (V4 - vi)k) • 

Here Mq is the Born amplitude with the extracted unstable propagators, Vi and V4 
stand for the top and the 6-quark three velocities respectively. As it is clear from 
this equation, the integration can be performed in a way similar to the case of the 
single resonance production. We introduce a proper time for each of the resonance 
propagators and exponentiate them. The integration over k is then the same as in 
the single resonance case (section 2). Finally we get: 

— «sCFl-ViVi /• (hdTi (./Di + D^ , Di 

IVln 



^part 



In what follows we denote the angle between velocity of the particle labelled i 
and the top quark velocity Vi by 9i. Then ri and r2 are: 

ri = 2/3 cos 6'4r — (1 — /3 cos ^4)ti 

and 

Next we make the change of the variables: t = Xx and Ti = X{l — x) with < a; < 1 
and < A < 00. Integrating further over A we obtain: 

ia.Cp , } dx 1 



We denote cos 9^ as X4 and : 



{I - j3xi) - x{l + I3xi) +ij?{l + x) 



i.'^%ni-x\). 

£/4 

Examining previous equations we see that the leading term under the square route 
can go through zero within the integration region if x^ > O.In this case this "would 
be" divergence is regularized by keeping the mass of the light particle finite. This 
means that the divergence is collinear. Actually this divergence can appear only 
if the "mass" of the gluon is zero, i.e. when the gluon pole in the original expres- 
sion is taken. Hence the appearance of this divergence in the particle pole means 
that decomposition of the unstable propagators (see eq.(34)) which we use for the 
evaluation of this graph has really mixed gluon and particle poles of the original 
expression in a nontrivial way. 
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We insert the identity 1 — 0{x4) + 9{—X4) inside the integral. After this we get: 

A = log{^l±^) + log(l - Px,) - log(-2/?a;4), 
B = log (1^^) - log [^^^) + log(l - (3x,) + log(2/3x4) 



log[/3^(l-a:/)]+log(^). 



^m4^ 

As it has been mentioned before and is quite clear from the above equation there 
are coUinear logarithms associated with the massless 6-quark in the final state. We 
discuss below (see section 3) how coUinear logarithms cancel in the final result. 

2.2.2 Gluon pole — Using decomposition eq.(34) for the resonance propagators 
and performing the integration over the contours discussed above we are forced to 
take the lower pole of the gluon propagator for the first term in the decomposition 
and the upper one for the second. Performing the integration over the modulus of 
the three-momentum we obtain the following representation for the amplitude: 

Where Ji and J2 are: 

J,^.[^^ ^(,og(H;^)+log(-^)-«). (37) 

J 27r4(-a;)(l - iifcns) ^ ^ 1 - (3x ^Di + D2^ ' 

Here is the unit vector, by x we denote cos 9k — n^ni and iij is the unit vector 
parallel to the velocity of the particle i. 

The integration over azimuthal angle is easily performed using the following 
equation: 

(r^= , (38) 



i l-n^n, y^(cos^fe-cos^i) 



This equation exhibits coUinear singularities which appear when the momentum 
of the gluon is parallel to the momentum of the (anti)quark. We regularize them 
keeping the mass of the light particle in the singular terms. The exact formula reads: 



\x - Xi\^ -^{x - Xi) + ^{l-Xi'^). (39) 

Finally changing the sign of the integration variable in the eq.(39) cosOk — >■ — cos^^ 
we get: 

Ji + J2=j f"" ( log (\^) + log (^) + ine{-x) - in9{x)) . (40) 
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It is rather straightforward to calculate the integral in the last equation. We 
split the integration region into two parts to rewrite the square root correctly and 
use a partial fractioning to obtain Spence-like integrals. 

It is quite useful here to examine a part of the previous expression which contains 
^-functions. The evaluation is straightforward. The result is: 



asCp (1 - V1V4) 
TT 2Di^{x4)m1 



A = {me{-Xi) - m9{xi)) Ai + 2m9{x^) A2 + 2i'Kd{-X4) A3. 
Where MgQ is the piece of the gluon pole part proportional to the ^-functions and 

--'o.@.l„,(fg)).l„,0), 



1 (Di + D2\ , /I +a;4\ 



^3 = -l0g(^^)-l0g(-^). 

If we sum the particle pole contribution elaborated above and the "-part 
of the contribution due to the gluon pole the result appears to be simple and all 
^-functions drop out. 

Mpart + Mge = Mo^-^^ !; , K, (42) 

TT Z U\C^{Xi) 

In this equation L4 stands for 

i4 = log(:^). 

This notation will be used further. 

To evaluate remaining contributions due to gluon pole it is convenient to use 
additional functions introduced in the Appendix. Finally we get the following result 
for the four-point function: 

agCp 1 — 0X4 . 
M,-, = Mo ' 2 ^tl. (43) 



27r Di^{x4)m\ 



A^l ^Ai+A2 + A3, 
Ai = -Fi{-Do, f3\x4) + Fi{-Do, -0\x4) + Fi{x4, Plx^) - Fi{x4, -f3\x4), 

A2 = log (^) \F2{x4x4) - F2{-Do\x4)\ , 
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In this equation we denote: 



It is also straightforward to consider gluon pole of the original matrix element 
(without decomposing resonance propagators eq.(34)). We need it due to the study 
of the bremsstrahlung integral, namely the interference of the gluon radiation from 
t and h quarks. The calculation is similar to the one described above. Finally we 
get: 

ITT OigCp 1 — 0X4 A , . .X 

Ag^Ai + A2, 

Ai = -Fi{-Do, f3\x4) + Fi{-Do, -/3|x4) + ^1(^4, P\x4) - Fi{x4, -P\x4), 

A2 = (log (^) + ZTt) {F2{x4\x4) - F2{-Do\x4)). 

Corresponding bremsstrahlung integral can be obtained from the previous equation 
by the standard change ( cf. discussion after eq.(32)) in the relevant piece of the 
differential cross-section. 



3.3 Five-point function 

We are finally left with the last non-factorizable graph which corresponds to the 
interaction between b and b. We again calculate a contribution due to particle and 
gluon poles separately. 

3.3.1. Pcirticle poles — Evaluation of the particle pole proceeds in a way similar 
to the one which has been used for the four-point function. We use decomposition for 
the propagators (cq.(34)) and then choose appropriate contour for the integration. 
The integral naturally splits into two pieces which represent the movement of a 
system of particles in the Coulomb field produced by a quark or an antiquark. 
These two pieces are symmetric and complementary to each other. 

Let us examine one of them. We exponentiate the propagators to obtain the 
Coulomb-like three momentum integral. As we have one more propagator here in 
comparison with the four-point function we need to introduce three "times" instead 
of two: 

Mj,artA = Mo 47rza,C^^i^^^ /, (45) 

/ = j drdndte^v ii-^r + -n) } J ^ (^^^^^ 

where r stands for the following vector: 

r = 2vT + (v - V4)ri + (vg - V4)t. (46) 

Here the quantity v is the on-shell velocity of the top quark. Integrating this equa- 
tion over k we get: 

. ^ (1-V3V4)^ f drdridt r,/Di + L>2 ^. 
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Here n4 is the unit vector parallel to the velocity of the particle 4 and is the 
component of the vector r perpendicular to the vector n4. 

The second term (Mpart.s) can be obtained from the eq.(52) by the following set 
of substitutions: 

V3 -V4, V4 -V3, Di D2, D2 Di. (48) 

The evident intention then is to perform the integration over t. The integral 
appears to be logarithmically divergent on the upper limit. This reflects the fact 
that infrared singularities of the five-point function can not be completely regularized 
by the virtualities and widths of the unstable particles. However, we anticipate that 
the above divergence corresponds to the Coulomb phase of the h quark in the field of 
the antiquark h. Hence we expect that this divergence is pure imaginary and drops 
from the observable quantities (as it occurs in the infrared limit of the "stable" 
theory |jl2|)- This expectation is verified by direct calculation. Below we omit this 
infinite piece from all expressions. 

The integration in the eq.(49) is then straightforward. We do not present its 
results because it is much more reasonable to present the sum of the particle pole 
and the -terms from the gluon pole. 

3.3.2 Gluon poles — Let us discuss the contribution due to gluon poles. As we 
have chosen appropriate contour to evaluate particle poles we are forced to take the 
lower and the upper poles in the gluon propagator for the first and the second term 
in the eq.(|3^) respectively. As in the case of the four-point function we perform the 
integration over the modulus of the three momentum and get: 

M, = -Mo^(l-n3n4) / ^^^f -^{D,, D^^cosd,). (49) 

2tx J (1 - nfen3)(l - nfen4) 

The function ^ can be written in the following way: 

^(Di, D2, x) = ^0(^1, ^2, x) + ^i(Di, D2, /3, x) + ^1(^2, D,, -f3, x) + (50) 



where: 



^^{D^, D2, P, x) = ^'^"^ , log ^ 



Di^{x)mf ^(1 — /3a;)m|^' 
2iTr(3x{e{x) - e{-x)) 



m1i{x){D^ + D2) ■ 

The parameter e in this equation is the infrared cut-off. The first term in the previous 
equation does not contribute to the observable quantities. Indeed, the infrared log 
is canceled by the real emission while the m term is pure imaginary and hence does 
not interfere with the Born amplitude. 

Next we evaluate the integral in the eq.(51). As the function ^ does not depend 
on the azimuthal angle, we calculate the following integral: 



'34 



. (51) 

27r(l - nfcn3)(l - ntn^) 
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Direct integration gives: 



(52) 



J ^ ^ 

(1 - n3n4)(l +X34)(x - Xa)(x - Xfe)' 

. N3-xKs N4-XK4 

34 — —I \ \ \ 

COS 9s + COS 64 ±i sin ^3 sin 64 sin (ps4 

^ ' 1 + cos 6'34 

A^3 = 1 — COS ^34 — COS ^3 (cos 63 — COS 64), 

N4 — 1 — COS ^34 — COS 6'4(cOS 64 — COS 6*3), 

K3 = COS 64 — COS COS ^^34, 

K4 = COS 6*3 — COS 6*4 COS ^34. 

Here we denote by a; = cos 6, by Xi = cos 6i and by X34 = cos ^34. Here the angle 
^34 is the angle between the vectors 03 and 114. 

As it is seen from this equation, the result of the azimuthal integration is a 
rational function of the cos(^). The remainder of the integrand consists of logs and 
constants, hence it is quite clear that the integration can be performed in terms of 
the Spence functions and logarithms. 

The other point is that divergence which occurs for x = 0:3,4 is the coUinear one 
and hence its regularization is clear. Explicit formula reads: 



\x 




We write eq.(54) in the following way: 



^34=:p— — (/3(X)+/4(X)) (53) 

i — n3n4 



where 

N. - xK. 



{1 + X34){x - Xa){x - Xb)\x - Xsl' ^^^^ 

Let us study the ^-terms of the ^ function and show how they cancel against cor- 
responding parts of the particle pole. As both particle and gluon poles are naturally 
splitted into two terms (3 and 4) we present them separately. Further evaluation is 
straightforward. The sum of the particle pole contribution and the ^-terms of the 
gluon pole reads: 

M41 = Mpart,4 + Mg^oA = -Mo^— K, (55) 

D2D ^irrif'' D2D1 ^imf^ 
2{l-f3x4), ^ N 2(1+/?X4), (D 



Di^{x4)m1 ^ ' D2^{x4)mf ^Di 
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i.{«42.og(%l)+.„g(|«))+ 

^i(log(l - ^i) - log(-l - Xi) - 21og(-a;i) + 21og(a;4 - + 

i=± 

R^(2 log(/3) + 21og(l-a;4)-L4)}. 



and 



M31 = Mpart,Z + My^efi 



— iasCp 



(56) 



2 / D 
log 



DiD 



D2D, 



log 



D2_ 



i{«3,(2.og(«f!i)..og@)- 

^i?i(log(l -Xj) -log(-l -Xj) + 21og(-a;i) - 21og(a;3 - 2;^)) + 



2(1 -/Sxs) 



+ 
D 



i?3(21og(/3) + 21og(l + X3)-L3)}. 

We denote D — Di + D2 in the above expression . The exact expressions for the 
quantities Ri can be found in the Appendix. 

Then we are left with the integration of the ^1 function. The result of the 
integration is: 
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-r-p OisCp 

-Mo 



E 



27r 

(1- 



Di^{x4)mi 



(log ( — y) (^) ~ ^1(^4, -P\x4 



(3xi 



(log (^^)F2(a;i|a;4) - Fi(a;i, -/3|a;4)) - 



TTlVaog (§)^2(-/^ok4) - F^i-Do, -(3\x,) 



+ 



{D^^D2,(3^-P). (57) 

Similar term (M32) which corresponds to the particle 3 can be then obtained by the 
direct substitution 3 4 in the eq.(59). Our final result for the radiative correction 
due to the bb interaction can be constructed from the above quantities: 



M, 



bb 



M41 + M42 + M31 + M32. 



(58) 



Finally we present the contribution of the "true" gluon pole (i.e. without de- 
composition of the resonance propagators eq.(34)) of the virtual five point function: 



Mg = M42 + M32 + M34 



(59) 



M34 = -Mo^-^ 



■iog(^)-E 

^ m4^ ^ ^ 



(1 ^/3X4) 

^-^F2(-Do|x4)] + (4^3). 
This concludes our evaluation of the five-point function. 



(1 - I3xi 



Ra 
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4 Analyses of the general formulae 



So far we have derived general formulae for the double resonance radiative correc- 
tions to the matrix element of the production of two resonances. Here we discuss 
some general properties of the obtained formulae. 



4.1 Collinear singularities 

As it is clearly seen from the above formulae each of the separate contributions to 
the non-factorizable radiative corrections exhibits collinear logarithms. Normally 
these logarithms are cancelled against the real emission. Let us note that T. D. Lee 
and M. Nauenberg [|T7| have used quite general approach to prove the absence of the 
similar divergencies in any quantum mechanical system. The basis for the proof is 
the existence of the unitary S-matrix. As is well-known from the work by M. Veltman 
||18| , it is indeed possible to construct the unitary S-matrix in the field theory with 
the unstable particle. Hence the arguments of the ref. Jl^ must apply also here. 



However it is necessary to clarify the level of the inclusiveness which is necessary for 
this cancellation to occur when unstable particles are considered. 

For this aim we extract all the terms which are singular in the limit — 0, 
2 = 3,4 from the above formulae and calculate their contribution to the total cross- 
section. In spite the fact that these terms are quite complicated in the individual 
graphs the sum of all these contributions appear to be very simple. We first write 
its contribution to the differential cross section: 

^ = 2 ^ Re| m (^L4 + ^La) I . (60) 

We remind that the quantities L3, L4 are defined by the eq.(42). 

Let us discuss now the properties of this equation. First we note that the source 
of this large logarithms are the virtual contribution due to the five-point function. 
Of course there are collinear logarimths also in the real interference but these are 
cancelled against similar pieces in the virtual corrections. 

We can also reexpress the terms in the eq. ( |60D to indicate exactly the mass 
singularities which we find in this case: 



dan 271 ^ ^ D D ^ ^ m 
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All other terms which have been dropped in the transition from the eq. (30) to 
the eq. (^) are smooth in the limit when the masses m^, i = 3, 4 go to zero. 

We see therefore that if the masses of light particles in the final state are equal 
then we do not get any mass singularities. This is the case for instance for the 
reaction e"'"e~ tt ^ W^W~hh. However such mass singularities will appear in 
the reactions like e~e^ W^W~ e^z/e/x+z/^. They drop out if the sum of charge 
conjugate channels is considered simultaneously - for instance e~e^ —>■ W~^W~ 
e~PefJ'~^J^fi and e~e'^ — > W~^W~ — > /i~z/^e"'"z/e. In any case if the integration over 
invariant masses of the produced resonances is performed, these singularities drop 
out from the observable cross section. 
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4.2 Real emission and virtual gluon poles 

We discuss here how the cancellation of the real emission and the virtual corrections 
occurs when unstable particles are produced. We begin with the single resonance 
production (cf. section 1). 

First we examine Be interaction (in terms of the section 1 ). It is straightforward 
to write the cross-section for the real emission and the contribution of the virtual 
photon pole to the cross-section in the soft photon approximation: 

«rt , remiss A^^\M \^A„^ f J ^ (Pl + k) + D jPl - k) ^ 

It is seen from this expression that gluon momentum enters the propagator of the 
unstable particle with different signs in virtual and real corrections. This is the illus- 
tration of the statement in the ref. 0] where the authors claim that the cancellation 
is not local in the momentum space in contrast to the usual situation. However 
the above expression is well-defined and we can evaluate it explicitly. The result of 
this calculation appears to be pure imaginary and hence does not contribute to the 
cross-section. The same situation also occurs for the usual triangle graph with one 
unstable particle. However the case with two resonances appears to be much more 
unusual. 

As is well known the usual thing in dealing with the soft limit of the Feynman 
graphs is the cancellation between real and virtual corrections. The essence of this 
cancellation is the fact that the particle movement is not affected by emission and 
absorbtion of soft massless quanta. Therefore the probability of a process remains 
the same. The piece of the virtual corrections that cancels real emission is the 
residue of the massless gauge boson propagator ( photon or gluon ). 

This simple remark verifies similar cancellation in the case when the integration 
over invariant masses of the unstable particles has been performed. In this case, as 
it is clear from our consideration, we effectively recover the situation with the stable 
particles. However the differential distributions represent a different case. 

Explicit investigation of the contribution due to the gluon pole from the virtual 
correction and the real emission shows that they cancel each other in a non-trivial 
way. Let us fix the off-shellness eq.(0) of one of the resonances Si. Then the virtual 
gluon pole contribution calculated for the off-shellness of the other resonance S2 
cancels the real emission for the off-shellness —62- The reason is that for negative 
values of 6 the particle is more likely to absorb gluons (the particle prefers to make 
its invariant mass larger) while for positive 5's the situation is opposite; exactly 
on the mass shell 6 = there is no difference. This shows that in the case of the 
unstable particle we have one more degree of freedom - the invariant mass which 
is sensitive to the soft {k ~ F) emission and absorbtion. Averaging over invariant 
masses we "lose" this degree of freedom (technically non-local cancellation in the 
space of the invariant masses occurs), but when the distribution in the invariant 
mass is studied we meet some unusual properties. 
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5 Conclusions 



We have derived general formulae for the non-factorizable radiative corrections to the 
invariant mass distributions for both single and double resonance production. We 
find these corrections to be important for the accurate description of this distribution 
in the vicinity of the resonance peak. 

Our approach is motivated by the observation that non-factorizable corrections 
to the Born amplitude are governed by the soft limit in order to give resonant 
contributions. This fact justifies the use of the soft photon (gluon) approximation 
for this problem. As usual the soft photon approximation provides universal results 
in the sense that they are not restricted to a concrete process. 

The gauge-invariant current can only be constructed if one takes into account 
both the current of the resonance and the current of its decay products. Gauge 
invariance is responsible for the cancellation of the whole effect for high energies and 
the most probable kinematical configuration, i.e. when the charge decay products 
follow the direction of motion of the resonance. 

We hope that our study provides better understanding of the structure of the 
infrared limit of the theories with the unstable particles. We note that the usual 
cancellation between soft real and virtual corrections is not complete even in the well 
known theories like QED: in fact the "photon" poles from the virtual corrections 
cancel the real emission, while the "particle" poles (which also give infrared diver- 
gencies) appear to be pure imaginary and physically correspond to the Coulomb 



phase H, [y 



In the case when we deal with the unstable particles the "particle" poles provide 
non-vanishing corrections to the observable quantities. The origin of this correction 
is very simple: the decay of the resonance accidentally changes the movement of 
the charge and hence destroys a coherence necessary to acquire "proper" Coulomb 
phase. Dealing with the Born amplitudes describing resonance production we can 
recognize that the integration over invariant masses of the resonance restores the 
"stable particle scenario" . As for the non-factorizable radiative corrections we know 
that they disappear if the integration over invariant masses is performed [4-6]. As the 
integration over invariant masses restores the stable particle scenario, the absence 
of the contribution due to non-factorizable corrections in the integrated quantities 
is in accordance with the non-observability of the Coulomb phase in the familiar 
theories with the stable particles. 

As for the cancellation of the virtual photon poles against the real emission we 
argue that this cancellation occurs only if the integration over invariant mass of at 
least one of the resonances is performed. 

Let us give a summary of the formulae presented in the text: 

• Non-factorizable radiative correction to the differential cross-section for a sin- 
gle resonance production is given by eq. fpl) . 



Non-factorizable corrections for the matrix elements describing production of 
two resonances are given by: 



1. three-point function - eq. 

2. four-point function - eq.(43) 
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3. five-point function - eq.(^). 

From the phenomenological side our study is motivated by a future investigation 
of heavy unstable particles. The vivid example is provided by the study of the 
reaction e~e~^ W~W^ at LEP 2. It seems that the planning accuracy of the 
measurement and the proposed technique of measuring the line shape of the invariant 
mass distribution requires taking into account QED non-factorizable corrections as 
well. 

As the dominant contribution to e~e^ W~W^ comes from the t-channel 
neutrino exchange, on the first glance it seems that the six-point function is actu- 
ally needed for this case. However, simple estimates show that practically for the 
whole phase-space of the final particles factorization of the Born amplitude is still 
valid. Hence it is sufficient to use the five-point function for the description of the 
production of two W bosons at LEP2. 

We note that the energy region for the LEP 2 {^/s = 170 — 200 GeV) is the inter- 
mediate but not really threshold energy region. Consequently one has to consider 
the effects of the final state interaction between decay products of the resonances as 
well: it is likely that the Coulomb correction alone (which is definitely the leading 
one in the threshold region) is not sufficient for the LEP2 energy region. 

As it has been indicated above our results beeing obtained in the soft approx- 
imation are universal. For the illustrative purposes we apply them to the process 
e'^e^ — » W^W~ — i> e'^i'e^'Ue- The values of different sources of virtual contributions 
as well as corresponding pieces in real interference are presented in the figs. 5-6. 

We stress however that these numerical consequences of our results for the ob- 
servable quantities seem to depend strongly on the experimental procedure which 
will be used in the real life experiments. It must be clear from the above discussion 
that the cancellation of the soft real emission against the virtual corrections is quite 
delicate in the case of the production of the unstable particles. Therefore a more 
realistic treatment of the soft {uj ~ r)radiation is necessary. This is basically the 
main reason why we do not see much sense in an exhaustive numerical analyses of 
our formulae. 

Another phenomenological issue which we mention here is the possibility to mea- 
sure the invariant mass distribution of the top quark at e^e~ and 77 colliders. Our 
formulae can be also applied for the 0(a<j) non-factorizable corrections in this case 
in the spirit of 0, |Tl|. Let us note however that in this case the infiuence of the 



hadronization on the precise determination of the top mass should be considered. 
The discussion of this important issue can be found in the ref. |T9[| . 

To conclude, we want to emphasize once more that the non-factorizable correc- 
tions change the shape of the invariant- mass distribution while preserve the total 
probability [4-6]. As the study of the properties of the unstable fundamental parti- 
cles requires the measurement of the invariant mass distributions our results must be 
taken into account while preparing for the analysis of the forthcoming high-precision 
experiments on the unstable particle production. Not only high statistics will be 
important but also our possibilities to make the correct correspondence between 
the results of the perturbative calculations of the masses and widths of the unstable 
particles in the framework of the Standard Model with the experimentally measured 
quantities. 
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Appendix 

Let us introduce the following integral: 

/dx 
—— log(l + (5x){e{x - Xi) - e{xi - x)). (62) 
X a 
-1 

Here a is a complex number with the non-zero imaginary part. Also Xi is an arbitrary 
number satisfying — 1 < < 1. The result of the integration then reads: 

Fi(a,/3|x,) = -21og(^^^|l^)log(l + /3x,) + log (ii±^) log(l - /3) + 

'os(^)Ml.«-2L,(i±$) . L,(i^^L,(i±|). 

In our formulas we also need this function in the case when a is real but in some 
restricted cases, namely a — Xi. In this special case the divergence is of the coUinear 
origin and we regularize it keeping the mass of the light particle finite. Hence the 
result for this function with a — Xi reads: 

F,(..,/3|.0 = log(l + /3x.)log(^) -Li,(^^-^) -Li,(^^-^). 

Our next function is defined as following: 

F2{a\xi) = [ -^{9{x - Xi) - 9{x, - x)). (63) 
J X — a 
-1 

The result of the integration is: 

F2{a\xi) = -2 log{xi -a) + log(-l - a) + log(l - a). (64) 
When a — Xi this function equals: 

F2ixi\xi)^log(^). (65) 
Next we present the quantities necessary for the eqs. (57-59): 

""^'^ mf(l + X34)e(^4)e(^+)e(^-) 

where all notations are the same as in the main text of the paper. The quantities 
for the index 3 can be obtained from the previous ones by direct substitution 4 — > 3. 
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Figure 1: Born graph and graphs responsible for the non-factorizable corrections for 
the simple model (see sect. 2). 
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Figure 2: The relative size of the non-factorizable radiative corrections in the simple 
model (see eq.(ll) with r] from eq.(15)). Curves A, B, C correspond to the total 
energies y/s — 180, 190, 200 GeV respectively. We use mw — 80 GeV and a — 
1/137. 



29 



0.04-- 



0.02-- 



-0.02-- 



-0.04-- 




Figure 5: Relative non-factorizable corrections to completely differential cross sec- 
tionon of the process e^e~ — > W^W^ — e^v^e^Vf. as a function of invariant mass of 
the e'^Ve system m2 in GeV for the fixed invariant mass of e~Pe "^i — 78 GeV. We use 
= 180 GeV, mw = 80 GeV, a = 1/137, Ow-e- = 30°, %-e+ = 150°, (^e+e- = 0. 
Curves A, B, C correspond to the conrtibutions due to three-, four-, and five-point 
functions respectively. 
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Figure 6: The same as in fig.4, but for mi = 82 GeV. 
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This figure "fig2-l.png" is available in "png" format from: 
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